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Exercice l (Discretisation l’Equation de transport (8 points)) 

On considere Vequation de transport lineaire a coefficient constant (c^O)en dimension 1 


du du 

Si +C ai =0 ’ (*.‘)SRXR+ 

u(x, o) = u°(x) , xen 

On approche Vequation parle schema suivant : 


u7+ ' ~ “ r ‘ + I ~ ” r ‘ + ^ ~ = 0 


2At 


(1) 


(2) 


On pose dans la suite a — 

h 

Question 1. Ce schema est-il implicite ou explicite ? 

Question 2. Calculer l’erreur de troncature de la methode. 

Question 3. Analyser la stabilite du schema (2) par la methode de Fourier. 

Question 4. Indiquer les avantages et les inconvenients quant a l’utilisation pratique de ce 
schema. 

On considere maintenant le schema numerique suivant : 


u ? +1 -^ , c . «? ' -»U\ _ n 

At + 2 1 h + h ) " ° 


(3) 


Question 5. Ce schema est-il implicite ou explicite ? 


Question 6. Montrer que le schema est consistant a l’ordre 2 sauf pour deux valeurs du pa 
rametres a (valeurs que Von precisera) pour lesquelles il est plus precis. 
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Question 7 . Analyser la stabilite du schema (3) par la methode de Fourier 

Question 8. Indiquer les avantages et les inconvenients quant a l’utilisation pratique de ce 
schema. Comparer au premier schema propose. 


Exercicb 2 (Modelisation du trafic routier (1 1 POINTS)) 

On modelise le trafic routier comme un flux continu de voitures par analogie avec les fluides 
compressibles. La densite de voitures (nombre de voitures par metre) est donnee par une fonc- 
tion u(x, t) dehnie sur R x R + . La vitesse instantanee est modelisee parv(x, t) — 1 - u(x, t). On 
suppose que u max = v max = 1. Ceci correspond a supposer que lorsque u = 0 la route est vide 
et les voitures ont une vitesse maximale. D’un autre cote , lorsque u = 1 ilyaun bouchon et les 
voitures sont a V arret. 


Pour une condition 
verihe : 


initiale de densite dehnie par u(x, 0) = u°(x), fonction de R — > [0, 1], u 



d(uv) 


dx 

u{x, 0) — u°(x). 


= 0, a: G R, i > 0 


(4) 


Question 1. Reecrire (4) sous sa forme conservative. Quelle est la fonction flux f(u) ? La fonc- 
tion f est-elle concave ou convexe ? 


Question 2. Montrer que la caracteristique qui passe par le point (£, 0) est donnee par 

*(0 = ( l-2u°(0)i + £ 

et montrer que la solution est constante le long de ces droites. 

Question 3 . Au travers d'une ligne de choc (u g , ua, £ = 7(i)), quelle est la relation de Rankine 
Hugoniot, liant la vitesse du choc aux valeurs de la solution a gauche (u g pour x < 7 (t)) et a 
droite (ud pourx > l(t))- 

Question 4 . Quelle est la condition pour ce probleme qui permet de determiner si un choc 
( u g , Ud , £ = 7 (t)) est entropique ? 


Question 5 . On considere la condition initiale : 


u°(x) = 


1 

- six < 0, 
1 six> 0. 


(5) 


(a) Tracer (de maniere approximative) les caracteristiques associees au probleme dans le 
demi-plan (x, t) e R x R + . Que se passe-t-il ? 

(b) Construire la solution faible entropique du probleme. Justiher votre reponse. 
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(c) On notep(t), la position d’une voiture a l’instant t 
dp(t) 


dt 


t), avecp(O) — £ la position initiale de cette voiture. 


En fonction des valeurs de determiner et tracer (de maniere approximative) la trajectoire p(t) 
dune voiture au cours du temps dans le demi-plan (x, t) e RxR+.En deduire une interpretation 
physique de 1 evolution du trade sous cette condition initiale. 

Question 6. On considere maintenant la condition initiale : 


u°(x) 


f . six < 0, 
1 

- si X > 0. 

O 


( 6 ) 


(a) Construire la solution faible entropique du probleme. Justifier votre reponse. 

(b) En fonction des valeurs de £, la position initiale d’une voiture, determiner et tracer (de 
maniere approximative) la trajectoire p(t) de cette voiture au cours du temps dans le demi-plan 

Or, t) 6 R x R .En deduire une interpretation physique de revolution du trade sous cette condi- 
tion initiale. 


Question 7 . On considere la condition initiale : 


u°(x ) - 


1 six < 1, 
1 

- Six > 1. 
z 


(7) 


(a) Tracer (de maniere approximative) les caracteristiques associees au probleme dans le 
demi-plan (x, t) e M x R + . Peut-on determiner u partout a l’aide de ces caracteristiques ? 

(b) Quelles sont les fonctions non constantes de la forme 


u(x, t) 




qui sont solutions au sens classique de (9) dans le demi-plan R x R+ ? 

(c) En deduire la construction d’une solution continue de (9) sur R x R+* associee a cette 
donnee initiale. Cette solution est-elle la solution faible entropique du probleme ? 

(d) En fonction des valeurs de £, la position initiale d’une voiture, determiner et tracer (de 
maniere approximative) la trajectoire p(t) de cette voiture au cours du temps dans le demi-plan 

(x, t) eRxR+.En deduire une interpretation physique de revolution du trafic sous cette condi- 
tion initiale. 


Question 8. On considere la condition initiale : 


u°(x) - < 


- six < 0, 

1 sj'O < x < 1, 


- six > 1. 


( 8 ) 
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Donnerla solution entropique du probleme. En deduire une interpretation physique de revolution 
du trafic sous cette condition initiale. 


Exercice 3 (Equations lin£arisees de la dynamique des gaz (7 points)) 

L ecoulement unidimensionnel d’un fluide parfait compressible est modelise, en coordonnees 
Euleriennes, par 

dp d 

+ ~dx'^ m ' = 0 , x 6 R, f > 0 (conservation de la masse) ( 9 ) 

d d 2 

Ot P u + +P) = 0, ieK,(>0 (conservation dela quantite demouvement) (10) 

9 d . . 

~Q t P e + fa. \\P e + P) u ) = 0 , x e R, t > 0 (conservation del’energie) (H) 

oil p represente la densite du fluide, u la vitesse, e l’energie totale specifique etp verifie la loi de 
comportement des gaz parfait pour une transformation isotherme 


p = Cp. 

Les conditions initiales sont donnees par 

p(x,0) = p°(x) u(x,0) = u°(x) e(x, 0) = e°(x) (12) 

L’etude de ce systeme de trois equations et trois inconnues presente de nombreuses difficultes 
theoriques dont certaines n’ontpas encore ete resolues a ce jour. Nous cherchons dans cet exer- 
cice a etudier ce probleme autour d’un etat d’equilibre. 


Question 1 On considere une solution (poo, Woo, eoo) stationnaire (independant du temps t) et 
uniforme (independant de la variable d’espacex ) du systeme (9) avec Poo > 0. On considere une 
perturbation autour de cet etat d’equilibre sous la forme 

(p(x, t),u(x,t),e(x,t)) = (poo^ocCoo) + e(pi(x, t),u\(x, t),e\(x,t)) + o(e) 

ecrire formellement le systeme verifie par (pi(x, t), ui(x, t), ei(x, t)) lorsquee tendversO. Onap- 
pelle ce systeme le linearise de (9) autour de (p^, Uoo , eoo ). 


Question 2. Ecrire ce systeme sous la forme d’un systeme du premier ordre 


d d 

jr.U + A—U = 0 oil A 
dt dx 


Est-ce que ce systeme est hyperbolique ? 


f u oo 

u oo 

0 \ 


' Pi(x,t)/p OQ ' 

\c/ Uoo 

'Moo 

0 

et U(x,t) = 

U\ (x, t )/ Uq q 

\ ^oo 

V'oo 'U'oc ) 


. e\(x, t)/C . 


(13) 
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Question 3. On pose 


Pi L Vc U\ . n , ^OO \ ei 

<p = 1 V (1 + — ==) — 

Poo V'OC V'OC VC C 

^ _ Pi VC ui | Uoo £1 

Poo ^OO ^00 \/C C 

Montrer que ces 2 fonctions ainsi que ei verident une equation de transport que Von precisera. 

Question 4. Calculer <t>(x, t), t ) et e\{x, t ) pour tout (x, t) e K x R + par la methode des 

caracteristiques en fonction dep\, w, ete\. En deduire p\ (x, t) etu\(x,t) pour tout (x,t) 6 RxR+. 

Question5. Trouver le plus petit intervalle[a(t), /3(t)] tel que si supp pj(-) c [-a, a], suppu ?(•) c 
[-a, a] et suppe^(-) C [ -a,a]alors suppp(-,t ) C [a(t),/3(t) ], suppu(-,t) C [a(t),/3(t)] et suppe^-, t) C 
[a(t),/3(Q] 
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